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04 ' We consider the application of the theory of vibrational control to _ffoo-problems. We study the possibility 

f— I ' of introduction of high-frequency parametric vibrations in order to decrease the minimal attainable value 

, of the iJoo-norm. We prove the existence of the stabilizing solution of the Riccati equation with quickly 

oscillating coefficients. This solution is found using the averaging technique as a series of the small parameter. 

One of the goals of control in various dynamical systems is their stabilization 
; with respect to all or to a part of phase variables. Traditionally, the main tool 

O ; of automatic control is feedback control, which is connected with the necessity 

to measure the state of the system. However, in some systems the measurement 
of the state is impossible. 

As an alternative to feedback control the principle of vibrational control was 
proposed by S. M. Meerkov M . The essence of vibrational control is in introduc- 



! tion of parametric vibrations into the system in order to change its dynamical 

O ■ properties. By means of vibrational control it is sometimes possible to achieve 



O 



the desired change without measuring the system's state. As an example of such 
^ ; system we can mention the inverted pendulum with the oscillating suspension 

^ ! point . The stabilization of the unstable equilibrium of the pendulum in this 

^ ' case occurs for sufficiently high frequency of vibrations. 

In the present article we consider the application of the vibrational control 
theory [1-4] to problems [5-8]. 
^ ! The general iifoo-problem in the linear system can be formulated as follows. 

Consider the linear system 

X = A{t)x + Bi{t)u + B2{t)w, 
z = Lxj 

x(0) = 0, tG[0,oo), 

where x G R", z G R"", u G Rp, w G R^^, A{t), Bi{t), B2{t) are continuous 
bounded on [0, 00) matrix functions. We assume that the pair {A{t), Bi{t)) is 
internally stabilizable, i.e. there exists a such continuous matrix function K{t) 
that the system 

x= [A{t) + Bi{t)K{t)]x 

is asymptotically stable. 



1 



Consider the space 

£^[0, oo) = {/ : [0, oo) ^ : /f (/(^), f{t))dt < oo 

and similarly defined spaces C2[^,oo), CHO^oo). In the sequel we shall denote 
these spaces as £2^, £2? ^l- 

Let U C £2 denote the set of controls of the form u = K{t)x internally 
stabihzing the pair {A{t), Bi{t)). We shall call the quantity 



2 ^Il7/I|2 
rm ~\~ \\Ui\\ rP 
L,2 II lli-2 



loo = sup — 2 



2 



the iifoo-norm of the system for the given u. 
Consider the functional 



J{u,w) = ll^;!!^^ + llwll^P - 7 11^11^.. (2) 
It is known that (e.g. 1^) that if the Riccati equation 

R = -A^R -RA + R [BiBI - ^BiB^^ R - LL^ (3) 
has the solution R^{t) = R^{t) > for t G [0, 00) such that the system 



X 



\ - [biBI - ^B2B^^ R, 



X 



is asymptotically stable, then 

1. The system x = [A — BiBiR^]x is asymptotically stable, 

2. The pair u* = —BiR^x, w* = ^B2R^x is such that 

J{u* ^w) = min J [u^w)^ J(w, w*) = max J(m, w;), J{u*^w*)=0. 

ueU weCl 

Moreover, if the solution R^ with the above mentioned properties exists, it 
constant if functions A, Bi, B2 are constant and is periodic if those functions 
are periodic. 

From 2 obviously follows that if such solution R^{t) exists, then J{u*, w) < 
for all w ^ jCl i-G- 



and 

l|-||oo<0. 

It is also known [0 that there exists 7* (J) - the smallest of 7 for which equation 
(^) has a solution. 



In the sequel we will show that 7* (J) can be decreased by the introduction 
of vibrations into the system. Assume that the system ([I|) has the form 

X = Ax + ^ sin ^Kx + Biu + B2W, , . 

z = Lx, 

where A, Bi, B2, K are constant matrices and e is a small parameter. 

As was mentioned before, the existence of a solution of the formulated above 
problem is related to the existence of positive definite periodic solution of the 
Riccati equation 

'^ = -(A + -sm-K]^ R-R(A + -sm-K]+ (5) 

dt \ £ £ J \ £ £ J 

■T 



R [BiBI - \b2BJ^ R - LL' 



T 

Next we prove that such a solution exists for small £ if there exists a positive 
definite solution of the averaged algebraic Riccati equation. 
We introduce notation 

D = BiBl - ^B^Bl, C = LL^. 

First, we consider the differential Riccati equation 
dR 

— = £ \-A^{t)R - RAit) + RD{t)R - Cit)] , (6) 
dt 

where A^ D, C are continuous T-periodic functions, D C are symmetric and 
£ is a small parameter. We shall look for a solution R^{t) = R^{t) > of the 
equation (g) such that the system 

x= [A{t) - D{t)R,{t)] 

is asymptotically stable. 

Consider the following spaces 

RS"^" - the space of symmetric real nxn matrices, 

CS^^" - the space of continuous symmetric T-periodic functions, 

CSZ^^" - subspace of CSJJ.^", consisting of functions with the zero average. 

Obviously, 

if the elements of RS"^" are considered as constant functions. 

Now we introduce the averaging operator Mq : CS^^" RS"^" defined by 
the equality 

MoF = F{t)dt. 



Clearly, Mq is a projector to the subspace RS"^". We also define the operator 
To on CSZJ^,^" 

TqF = F - MqF. 

For brevity we will denote MqF as F. 
Consider the algebraic Riccati equation 

- A^R -RA + RDR -(7 = 0. (7) 

Theorem 1 // equation (Qj has such solution Rq = Rq > that the matrix 
A — DRq is Hurwitz, than there exists such > that for any £ < equation 
(Qj has such a T -periodic symmetric positive definite solution R^{t), that the 
system 

X = [A{t) - D{t)R^{t)]x 

is asymptotically stable. 

This solution can he approximately found as 

R{N){t) = Ro + e{Ri + Ui{t)) + . . . + e^'iRN + UN{t)) + e^'+'UN+iit), 

where Rq, Rn G RS"^", Ui{t), . . . ,UN+i{t) G CSZ^^", and for some C > 
independent of e 

sup \\R,{t) - RiN)it)\\ < C£^+\ (8) 

Proof. Formally substituting R(^n) into Eq. (0), we have 

Elii + e^U2 + . . . + e'^+'tlN+i = 

e[-A^{t) {Ro + £{Ri + Ui{t)) + ...)- 
{Ro + £{Ri + Ui{t)) + ...)A{t)+ 
{Ro + £{Ri + Ui{t)) + . . .) D{t) {Ro + £{Ri + Ui{t)) + ...)- C{t)] . (9) 
Equating in (Q) the terms of order we have 
r/TT 

— i = -A^{t)Ro - RoA{t) + RoD{t)Ro - C{t). 
at 

Averaging the latter equation yields 

- A^Rq - RqA + RqDRo -C = (10) 



A^{s)Ro - RoA{s) + RoD{s)Ro - C{s)\ds 



According to the assumptions of the theorem, equation (|TD|) has a positive 
definite solution Rq G RS"^", for which the matrix A — DRq is Hurwitz. 



Equating the terms of order we have 

' = -A'{t){Ri + ni(t)) - + 



dt 

{Ri + I{i{t))D{t)RQ + RoD{t){Ri + Ui{t)). 
Averaging of the latter equality yields 

{A^ - DRofRi + Ri{A^ - DRo)+ (11) 

Mo[-A^{t)Ui{t) - Ui{t)A{t) + Ui{t)D{t)Ro + RoD{t)Ui{t)] = 0. 

Because the matrix A — DRq is Hurwitz, the Lyapunov equation ([TT|) has the 
unique solution Ri. 

The matrix 112 (t) is defined by 

Mt) = ll To [A'{s){Ri + ni(s)) - {Ri + Ui{s))A{s)+ 

{Ri + Ui{s))D{s)Ro + RoD{s){Ri + ni(s))] ds. 

Similarly we can find i?2, • • • , Rn and n3(i), . . . ,IlN+i{t). 
Next we prove the existence of solution satisfying the conditions of the 

present theorem. Obviously, 



dR 



dt 



-A' {t)R^N) - RiN)Ait) + R^N)D{t)R^N) - C{t)+ 



where g{t, e) - T is a periodic in t and uniformly bounded in e G (0, Eq] function. 
First we prove that the equation 

dp) 

— = £ [-(^ - DR^^^fe - e{A - DR^^)) + eve - E''+'g{t, e)] (12) 

has a T-periodic solution. Note that the function H{Q) = QDQ is such that 
for any a > there exists ^ > such that if ||0i|| < S and ||0i|| < S, than 

||//(ei- 62)11 <a||ei-e2||. 

Let r = MqQ and A = TqQ. Than applying operators Mq and Tq to the left 
and right hand sides of Eq. ( [T^ ) , we have 

= [-{A - DR^N)f{r + A) - (r + A)(^ - Di?(A.))+ (13) 

{A - DR^N))D{A - DR^N)) - e^'+'git, e)] dt, 
A = eTo [-{A - DR(^j^^f{r + A) - (E + A)(^ - DR^n))+ (14) 
{A - DR^N^)D{A - Di?(^)) - £^+i^(t, e) . 



From Eq. ( p^ ) we have 



Ar + BA + Gi(r,A,e) = 0, 



where the operators A : RS 
defined by the equahty 



nxn 



RS"^" and B : CSZ^^" 



RS' 



(15) 



are 



Ar = -{A - DRo)^r - r{A - dRq), 

BA = -Mo [{A{t) - D{t)RofA + A{A{t) - D{t)Ro) 



and 



Gi(r, A, e) = Mo [eDpi{t, e){r + A)+ 
£(r + A)Dpi{t, £) + (r + A)D(r + A) - e) 

where pi{t,£) = iR(N)it) - Ro)/£- 

Obviously, for any di > there exist £o and 5 such that if HFiUrs < ^? 
||r2||Rs < <5, ||Ai||csz < ^, IIA2IICSZ < S and £ < £0 than 

||Gi(ri, Ai,£) - Gi(r2, A2,£)||csz < ^(||ri - r2||Rs + ||Ai - A2||csz). (16) 

Since the matrix A — DRo is Hurwitz, the operator A is invertible and from 
([T5|) follows that 

r = -A-^BA - A-^Gi{T, A, e). (17) 
From Eq. (|T^) we have 

A = £G2(r,A,£), (18) 

where 

G'2(r,A,£) = 

fjo [-{A{s) - D{s)R^^){s)f{T + A(s))- 

{T + A{s)){A{s) - D{s)R^^^{s))+ 

(r + A(s))D(r + A(s)) - £^+V(s, £)] ds. 

(^2(1^5 ^5 s) satisfy a condition similar to (|TBD, namely that for any cr > there 
exist such £0 and 5 that if HFiUrs < 5, ||F2||rs < ||Ai||csz < ||A2||csz < ^ 
and £ < £0 than 

||G2(Fi, Ai,£) - G2(F2, A2,£)||csz < - F2||rs + ||Ai - A2||csz). (19) 

System (P7|), ([T8D can be rewritten as 

-A-^B 



F 
A 











T\ -A-iGi(F,A,£] 
Aj^l £G2(F,A,£) 



(20) 



or 



(F,A) = (^(F,A), 



where 



I [A I \ eG2{T,A,e 

The map cp acts in the space 

X = RS"^" X CSZJ^^", 

where we introduce the vector norm 

||(r,A)||x = (||r||Rs,||A||cszf. 

Next we show that there exists a baU in X, within which the map (p satisfies 
the generalized contraction principle JTTI, i.e. there exists the matrix K G R^^^ 
with non-negative elements and the spectral radius p{K) < 1 such that 

||^(ri, Ai) - ^(r2, A2)||x < i^ll(ri, Ai) - (r2, A2)||x, 

where the relation of inequality is understood in the following sense: 
it is said that x < y {x^y G R^) if Xi < yi for z = 1, 2. 
From ( |2D| ) and estimates ([TBD, ([T^) follows that if HFiUks < ^, ||r2||RS < 



|Ai||csz < (5, IIA2IICSZ < 5, then 

||^(ri,Ai)-^(r2,A2)||x< (21) 

£(72 £(72 J 

where di, (12 ^ for (5 ^ 0. 

Therefore, the map (p satisfies the generalized contraction principle for suffi- 
ciently small 6. Let us now find the value of S for which (p maps the ball of the 
radius 6 in the space X into itself. 

From (|T^) and (|r9|) follows that 



||Gi(r,A,£)-Gi(0,0,£)||<ai(||r|| + ||A 

||G2(r,A,£)-G2(o,o,£)|| <a2(||r|| + ||A 

for ||r|| < 6 ||A|| < 6. This implies that 

||Gi(r,A,£)||<Ci£^+i + ai(||r|| + ||A 
||G'2(r,A,£)||<C2£^+i + a2(||r|| + ||A 

where Ci, C2 are some constants. 



From (^) follows that 

ll^(r,A)||x< 



i.e. for some C* (f maps the ball of the radius p = C*e^^^ in the space X 
into itself and is contracting, and hence has in this ball the unique fixed point 
(r^, A^(t)) such that 6*(t) = + A^(t) is the solution of Eq. ([T2|). It is easy 
to verify that R^{t) = R(^N){t) + 6(t) is a solution of equation(|]). 
Next we prove the asymptotic stability of the system 

X = [A{t) - D{t)R4t)]x. (22) 

Substituting in the form Rq + £pi{t) + Q{t) into ( |22| ) and taking into 

account that the matrix A—DRq is Hurwitz, we have that (^) is asymptotically 
stable for small e. 

Positive definiteness of R^{t) for small e follows from the positive definiteness 
of Rq. 

Theorem 1 is proved. 

Consider Eq. (^. Introduce the fast time r = t/e. Then 

^ = - sin t(K^R + RK) + ei-A^R - RA + RDR - C). (23) 
dr 

Putting £ = 0, we have 

dR . , ^^j. 



sinriK^R + RK). (24) 
dr 

Equation ( p^l) has the general solution 

R = m{T)P'^^{T), (25) 

where P is an arbitrary constant, ^{t) = exp(A'^(cosr — 1)). 

Introducing P, defined by ([25]) , as the new variable in (^), we have 

dP . 



, =£[-A\t)P-PA{t) + PV{t)P-C), (26) 
dr ^ ' 

where 

A{t) = ^^{r)A{^^)-\T) V{r) = ^^(r)OT(r), 
CiT) = ^-\T)Ci^^)-\T). 
Consider the algebraic Riccati equation 

- A^P -PA + PVP - C = 0, (27) 

which is the averaging of the right hand side of 



Q 



Theorem 2 // equation ( ^ ) has a solution Pq = Pq > such that the matrix 
A — VPq is Hurwitz, then there exists a > such that for any £ < £q equation 
has a 27^ £ -periodic symmetric positive definite solution for which the 

zero solution of the system 

dx 



dt 



A' 



X 



is exponentially stable. 

For any N > this solution can be found in the form 



R 



{N) 



.£, 



.£, 



Pn + £ hPl + Hi - + • • • + 



£ 



N 



Pn + Un - +£ 



.£, 



where Pi, 



.£, 

rnxn 



Pn G RS"^", Hi, ... , n^+i G CSZ^^^", and for some C>0 

ft' 



sup 

Te[0,Te] 



R 



{N) 



.£, 



< C£ 



N+1 



Proof. According to Theorem 1, equation ([26|) has a positive definite stabi- 
fizing solution P^{t) such that for 



P{N){r) = Po + £{Pi + ni(r)) + . . . + £^(Pa. + UN{r)) + £ 



n 



■j) 



the estimate (g) is satisfied. It is obvious that R{t/£) = {t / £) P^{t / £)'^^ {t / £) 
is a solution of (|^) and for P(7v)(^/^) R*{t/£) a similar estimate is satisfied, 
and that the system 

x= [A{t) -V{r)R,{r)]x 

is exponentially stable. 

Example. Let matrices A, B, C and K have the form 



A 





-0.27 



1 

-2.8 



B 




1 



c = 



1 
1 



K 





k 



The averaged Riccati equation ([27| ) in this case has the form 



-0,27-A;V2 

1 -2,8 



P + P 



1 

-0,27-/c72 -2,8 




1 



p + 



1 
1 



= 0. 



+ 



n 



The table below shows the values 7^ for different k. These values are found nu- 
merically using the conditions for the existence of solutions of Riccati equations, 



given m 



k 







3,704 


0,25 


3,320 


0,5 


2,532 


0,75 


1,815 


1 


1,300 


1,25 


0,925 


1,5 


0,717 


1,75 


0,556 



This example shows the possibility of introduction of vibrations into a linear 
system, which decrease 7^ i.e. decrease the influence of disturbances on the 
output of the system. 
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